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FOREWORD 



This report describes work accomplished in the program "Quasi 
Area Rule for Heat Addition in Transonic and Supersonic Flight Regimes" 
conducted under USAF MIPR APO-71-007. The work was accomplished 
during the period 1 Oct 1970 through 31 Aug 1971 at the Naval Post- 
graduate School, Monterey, California, under the direction of Dr. Allen 
E. Fuhs. The report was submitted on 1 Nov 1971. 

The program was sponsored by the Air Force Aero Propulsion Labora- 
tory, Wright-Patterson Air Force Base, Ohio under Project 3066, Tur- 
bine Engine Propulsion, Task 306603, Advanced Component Research. 

Dr. Kervyn D. Mach, AFAPL/TBC, Turbine Engine Components Branch, 
was the project engineer. 

Publication of this report does not constitute Air Force ap- 
proval of the report's findings or conclusions. It is published 
only for the exchange and stimulation of ideas. 




ERNEST C. SIMPSON 

Director, Turbine Engine Division 



ii 



ABSTRACT 



Body shapes, including axisynmetric and three dimensional, have been 
developed to minimize wave drag. The von Karman ogive and the area rule are 
examples. Similar work has not been accomplished for optimum shapes with 
propulsion. Prop\ilsion can be divided into two categories — those devices with 
internal heat addition and those with external burning. For internal heat 
addition an analytical model is formulated which introduces the propulsive disc. 
Attention is shifted to external burning, which is examined for one dimensional 
and two dimensional linearized flow. Heat fronts and combustion fans are 
discussed as examples. Forces on a heat source in a uniform stream and 
adjacent to bodies are derived. Several possible applications are exjunined 
including base pressure augmentation by external burning, spin recovery using 
external burning, and transonic boattail drag alleviation. Previous work on 
base pressure augmentation has used a two dimensional planar model. A two 
dimensional eixisymmetric model is examined. 
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QUASI AREA RULE 



FOR HEAT ADDITION IN TRANSONIC 
AND SUPERSONIC FLIGHT REGIMES 

Dr. Allen E. Fuhs 
Professor of Aeronautics 
Naval Postgraduate School 

INTRODUCTION 

In subsonic flow, drag can be attributed to viscous forces on the body 
surface and to the release of vorticity in the wake. The vorticity is associated 
with the generation of lift. In supersonic flow, there is a new phenomenon lead- 
ing to drag; it is the radiation of waves from the body. These waves carry 
energy whose origin is the work done by drag and vehicle motion. 

For almost every case of a body flying through a medivtm, drag is considered 
to be a degradation of performance. Consequently, considerable effort has been 
devoted to optimum shapes giving minimum drag. 

One approach for optimizing a body shape is to represent the body by sources 
and sinks. For a closed axisymraetric body, all sources and sinks are on the 
axis. It is found that for a nonlifting body, the slope of the cross sectional 
area curve equals the distribution of sources and sinks, f(x). Distance along 
the axis is x. By means of the calculus of variations, one can arrive at shapes 
giving minimum wave drag. There are constraints on the body geometry, and the 
opt im u m shape depends on these constraints. Two examples are given in Table I. 
The shapes given in Table I are for minimum wave drag. Viscous effects are not 
included, and the bodies are nonlifting. 
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Table I. Optimim Shapes 



Name of Body Constraints Wave Drag 

Von Karman* Ogive Given length, t, and base area, 3 

Sears^^Ogive^^ ^ Given length, t, and volxime, V 

^‘Used on Low Altitude Short Range Missile. 

^'^Discovered by W. R. Sears of Cornell University. 




It was discovered by Lomax and Heaslet^^^ and by Whitcomb^ that the 
optimum shape for an axisymmetric body applies to slender bodies with wings and 
rudders if one interprets the area distribution properly. It is the area cut by 
the Mach cone. If one uses this area distribution from a three-dimensional body 
to match the area distribution of an optimum axisymmetric body, then minimum wave 
drag occurs. This procedure is known as the area rule. 

For a body with thrust, there are several new variables introduced. Thrust 
may be developed by adding heat to a flow internal to the body. Thrust also may 
be obtained by adding heat external to the body. 

Consider external burning. Is there an area rule for a body with external 
heat addition analogous to the area rule involving body shape only? For this case, 
there are two functions to optimize. One is body shape, and the other is the 
distribution of heat addition. Table I gives some typical constraints; an 
additional constraint for the optimization could be specific fuel consumption. 
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Instead of finding minimum drag, the procedure should give the best thrust- 
minus-drag or zero thrust-minus-drag depending on the constraints. Since both the 
distribution of heat and cross sectional area are involved, an appropriate name 
is "quasi area rule." 

When thrust is developed by adding heat to a flow internal to the body, the 
quasi area rule suggested in the previous paragraphs does not apply. For in- 
ternal burning, the waves due to heat addition obviously do not occur as in the 
case of external burning. There is, however, a change in stagnation pressure 
and temperature between the inlet and nozzle. Also the velocity vector changes 
from inlet to exhaust. An energy disc can represent those changes. This 
report discusses aspects of the thrust-minus-drag problem for an airbreathing 
engine with confined heat addition. Figure 1 illustrates the various problems 
that can be studied with internal burning when the emphasis is on thrust-minus- 
drag. The analysis is linearized. Two-dimensional (2D) planar is the geometry 
of a propulsive wing. Axisymmetric 2D geometry is an isolated nacelle without 
fins or wings. Addition of fins or wings yields the 2D propulsion with 3D body. 

The connection between these two geometries is the area rule. 

Figure 2 outlines the problems of interest for external burning. The 
analysis is linearized. For the 2D planar case, the heat and mass sources can 
be distributed over a surface or throughout a volume. Volume distributions must 
meet requirements for slender body theory. Axisymmetric 2D bodies are related 
to 3D bodies by the quasi area rule. 

The linearized 2D planar flow with heat addition has been studied extensively. 
A solution to the 2D axisymmetric case does not appear in the scientific 
literatxire. Solutions for 3D bodies with heat addition have not, as yet, been 
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Figure 2. Classification of Flows with External Burning. 









obtained. In order to develop the quasi area rule, these solutions are 
essential building blocks. 



PROPULSION WITH INTERNAL HEAT ADDITION 



Introduction 

For a supersonic aircraft, the variation in nozzle pressure ratio (NPR) 
is very large. This requires a variable geometry nozzle. Under certain flight 
conditions, the nozzle setting and nozzle flow result in large base or boattail 
drag. 

One motivation for treating the propulsion-airframe problem as a unit from 
the area rule point of view is to obtain new design approaches. First, the 
typical variation in NPR and nozzle geometry will be calculated for a 
Mach 2.5 aircraft. These results will be discussed to illustrate the problem. 
Second, the model for representation of a propulsion system with internal heat 
addition will be developed. Third, approaches to solution of the model will be 
outlined. 

Exhaust Variation with Mach Number 

The aim of this section is to calctilate nozzle pressure ratio, nozzle area 
ratio, mass flow ratio, and thrust of an ideal turbojet as a function of Mach 
number. Knowledge of the trends is essential in understanding aircraft exhaust 
problems. 

Consider an ideal turbojet with an ideal inlet. There is, however, a 
restriction on the Mach number, M^f into the compressor. The thrust is given by 



where m is a mass flow of air, and m, is mass flow of inlet at design conditions, 
a ’ d 




( 1 ) 
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For the inlet, the various areas are illustrated in Figure 3(a). Station 0 
is at a point in the freestream where is ambient pressure. Stations 1 and 2 
are inlet throat and compressor face respectively. The mass flow ratio is 



m 



m. 



^0 ^ \ 

A| 



( 2 ) 



“d ‘"od ^ H 

where Mq is flight Mach number, and M^^ is design flight Mach number. Areas 
Aq» A^, and A^^ are identified in Figure 3(a). Both a fixed area and variable 
geometry inlet are considered. The ratio A^/A^ is obtained from isentropic 
flow tables for given Mq for the variable geometry inlet. For a fixed geometry 
inlet Aq/A^ comes from flow tables, entering the flow tables with the Mach 
number downstream of a normal shock at Mq. The mass flow ratios are shown in 
Figure 4 for both cases. A variable geometry operates supercritically and 
swallows more air. A variable geometry ramp or spike inlet may have superior 
mass hgmdling capability at subsonic speeds. 

The thrust per pound of air for an ideal turbojet is given by 
F 



w 




T* 1 - 



T'r^T'r^i 
r ci 



- Tjr - 1) 

r c 



- 1 



(3) 



where the symbols not previously defined have the following meaning: 
a^ = freestream speed of sound 

7^ = T^q/Tq = ratio of stagnation to static freestream temperature 
7^ = T^^/T ^2 “ stagnation temperature ratio across the compressor 
The symbol T* deserves more discussion. It is 

T* ' ' V^O 



(4) 



where = stagnation temperature ratio across the burner or combustor. 

The value of T^^ is limited and has a maximiun value. Hence for a given engine 

and altitude, T* is fixed when the engine is at maximum power or maximum TIT, 
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FIGURE 4. MASS FLOW OF AIR AS A FUNCTION OF FLIGHT MACH NUMBER. 
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The station numbering is shovm in Figure ?(b). The throat area for the 



nozzle is 






V^’^TS "'a 



( 5 ) 



r Px8 

neglecting fuel addition. R is the gas constant. T is a function of the ratio of 
heat capacities, Y> has a value of ,683 for y ~ 1»4. The nozzle pressure ratio, 
NPR, is p^g/pQ. A ratio of Aj to the design value is 



A, 



7d NPR y T.^g^ 



(6) 



The area Ag is obtained from NPR, isentropic flow tables, and the assumption 

Pg = Pq. Two useful ratios are Ag/A^ and Ay/A^, i. e., comparing nozzle areas with 

compressor face area. 

Calculations were made for the values shown in Table II, 

One of the results of the calculation is NPR which is shown in Figure 5. Some 
of the other results are summarized in Table III The symbols for the turbine ratios 

are TJ. = '^t 5^^T4 "^t " Pt5^^T4* 

Figure 6 illustrates the area variations required for a supersonic turbojet 
operating at maximvim TIT for different Kq. Examination of Figure 6 and the data in 
Table III shows why turbojet nozzles usually operate underexpanded at design Mach 
nviraber. The percentage variation in Ay is 38 per cent based on design Ay. The 
percentage variation in Ag is 77 per cent based on design Ag. NPR varies from 4 at 
low subsonic flight to more than 25 at design of 2.5. 
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Table II. Values for Calculation of Turbojet Quantities 



Maximum turbine inlet temperature 
Altitude 

Freestream density 
Freestream speed of sound 
Compressor area 
Tip radius at face 
Hub radius at face 
Maximum Mach number at compressor 
Inlet throat area at design 
Design flight Mach number 
Compressor pressure ratio at 
100 per cent corrected speed 




PQ 

a« 



M, 



2max 



‘td 



M, 



Od 



2000°F 
Sea Level 
.002378 slugs/ft^ 
1117 ft/sec 

8.05 ft^ 

2 ft 
1.2 ft 
0.65 
7.1 ft^ 

2.5 



TT 8 

C 
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Table III. Results of Calculation for Turbojet Performance and Areas 
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As shown in Figure 6, to obtain ideal expansion the nozzle is converging- 



diverging for all Mq illustrated. Usually the nozzle is a simple converging 
nozzle from subsonic through transonic range of Mq, This results in a large 
boattail area. At the design Kq the exit area for ideal expansion significantly 
exceeds the maximum area of the propulsion system. The nozzle is cut off, relative 
to ideal, at some area appropriate to maximum engine area. The result is under- 
expansion and a ballooning plume, A shock wave may originate at the diverging 
streamlines of the plume, giving a drag increment. 

Model for Propulsion with Internal Heat Addition 

Consider a propulsion device as sketched in the top of Figure 7, To 
represent the wave drag of the inlet, it is necessary to duplicate streamline 
A-A, It does not matter what the details are of the inlet geometry internal 
to the streamtube bounded by A-A. 

The inner body which generates a flow having streamline A-A is shown. Note 
that it is not the same shape as the inlet spike, A distribution of sources on 
the axis could generate the inner body as well as streamline A-A. Note that the 
spike shock wave does not have the same angle as the cowl shock wave. 

If one applies the momentum relations to the control volume of Figure 7, 
the thrust minus drag is given by 
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T - D 1 



Dr 

(p„ - p )dA - — + 2 
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FIGURE 7. MODEL OF AN INLET USING SOURCES. 



where and are the dimensionless perturbation velocities arising 

from the nacelle. It is apparent that for the case where the exhaust exit 
defines a plane, as illustrated in Figure the thrust and drag can be 
clearly separated. Waves generated by the exhaust plume do not influence T - D. 

Following the analysis suggested by the model of Figure 7, one can obtain 
a nacelle of known wave drag. Figure 8 illustrates the concept. Start with an 
ogive of known shape. Plot the streamlines in the flow adjacent to the ogive. 
Pick the radii rg and rg suitable for the nacelle. Insert the nacelle within 
the streamtube. The wave drag of the ogive and nacelle will be the same. The 
drag can be calculated using the potential function for the ogive. In the 
lower part of the figure, a nacelle with a spike inlet is shown. 

Consider some thrust device inside a momentum control volume shown in 
Figure 9. How does one modify the variables so as to change conditions on the 
surface of control volume and improve T - D? What does one have available to 
change the wave drag? Available variables are: (a) number and location of 

nozzles, (b) velocity profile at nozzle exit, and (c) body geometry. 

To obtain an optimum configuration, one must have a method to represent 
the body and exhaust analytically. The analytical representation can be 
subjected to the calculus of variations. 

The body can be represented by sources and sinks. Changes in stagnation 
pressure, stagnation temperature, velocity vector, and other flow quantities 
make representation of exhaust more difficult. 

Carrying the idea of Figure 7 further, let us look at an ideal ramjet. 

For an ideal ramjet, the following equalities hold true: 
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For an ideal ramjet (or for any other propulsion device satisfying the above set 
of equalities), one can represent the flow by sources, sinks, and an "energy 
disc." This is illustrated in Figure 10; stations 1 and 6 are identified in 
Figure 10. The inner body duplicates streamline ABC. The inner body is formed 
from the source-sink distribution in the right-hand column of Figure 10. The 
change in density, velocity, and stagnation temperature is achieved by an 
"energy disc." In fact, that defines the energy disc. It changes the flow so that 
the equalities list above is correct. Tne sources, sinks, and energy disc shown 
in the bottom of the right column would have the same thrust and drag as the 
ideal ramjet sketched in the left-hand column. 

This model looks as if it is satisfactory for an ideal ramjet, but it has 
some drawbacks: 

1. The energy disc is not a heat front or other simple fluid flow. 

2. When / A^, it may be difficult to apply the energy disc. 

1. For the case p^^^ / p^^^, some means of changing p^ must be added. 

The energy disc is bounded by a streamtube. Changing p^ within the disc gives a 
mismatch across the disc in regard to streamline shape. 

A satisfactory way of specifying the conditions at station 6 is needed, and 
the energy disc has been suggested. One wants to retain sources and sinks for 
two reasons. Wave drag is readily calculated if you know the distribution function.. 
The distribution function can be optimized to yield minimum wave drag. 
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FIGURE 10. IDEAL RAMJET REPRESENTED BY SOURCES, SINKS, AND AN ENERGY DISC. 



The energy disc represents the propulsion device. A particular engine 
generates changes in p^, T^, and velocity between the inlet and exhaust. For 
the case of a plane for an exit as shown at station 6 in Figure lOy one does not 
need the complexity of an energy disc. Equation (^7) indicates that drag and 
thrust are separable for this case. When one varies quantities (a), (b), and 
(c) discussed previously as available, independent variables, then the energy 
disc may be an attractive model. 

This section is closed by noting that a model has been proposed but 
has not been developed. 

EXTERNAL BURNING; BASIC EQUATIONS AND SOLUTIONS 
One-Dimensional Heat Addition 

To provide a basis for comparison with subsequent sections, one-dimensional 
heat addition will be briefly reviewed. 

There are several approaches to one-dimensional heat addition. For unsteady 

( 6 ) 

flow there is the self similar technique developed by Sedov. We consider 

steady flow only here. There is the famous solution for constant area heat 

addition, i.e., the solution by Rayleigh. Rayleigh flow is discussed in most 

( 7 ) 

gas dynamics books. When there is simultaneous variation of area and heat 
addition, the techniques outlined in the Princeton Series, Volume III, are 
applicable. 

( 9 ) 

The influence coefficients develoj^ed by Shapiro and Hawthorne will be 
discussed here. Table IV summarizes the influence coefficients for heat 
addition and isentropic flow with area change. The influence functions relate a 
fractional change of a dependent variable, say M^, to a fractional change of the 
independent variable, A or h^. So far as heat addition is concerned, the 
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Table IV, Influence 


Coefficients for Area Change and 


Heat Addition 




dA 


dh 




T 


~cl 






P 


dM^ 


-2(1 + 


1 + 




1 - 


2 

1 - M 


du 


1 


1 


u 


1 - 


1 - 


dT 


(y - 1)K^ 


1 - yM^ 


T 


1 - 


1 - 






1 


P 


1 - 


1 - 


d£ 




yM^ 


P 


1 - 


1 - 


dF 


1 






2 


0 


F 


1 ♦ ym 




ds 








0 


1 


P 







velocity, u, and static pressure, p, change as they would flow into decreas- 
ing area. For u and p, heat addition and area variation produce the same 
result. 

For itnconfined flow with heat addition, there are some analogies between 
the flow due to a particular body shape and the flow due to a given distribution 
of heat addition. These analogies, which are limited in scope, will be pointed 
out as they arise in subsequent discussion. The fact that heat addition and a 
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body may influence the flow the same way may be understood using the du/u and dp/p 
influence coefficients. 

The influence coefficients are the first terms in a Taylor series expansion 
of the dependent variables in terms of the independent variables. Many of the 
results of linearized heat addition in an unconfined flow can be related to 
influence coefficients. This statement is applicable, for examble, to heat 
fronts and combustion fans. 

Linearized Equations for Combined Mass and Heat Addition 

The equations of motion for steady flow with heat and mass addition are as 
follows : 



Continuity 


div(pu ) = m 


(8) 


Momentum 


Du ^ 

p^ VP + mu = 0 


(9) 




2 




Energy 


div piT(h ■•■ ^) = Q 


(10) 



In the preceding equations, m is the rate at which mass is added per unit 
volume. Likewise Q is the rate at which heat is added per unit volume. 

To linearize these equations, one introduces perturbation quantities which 
are defined as follows: 



p = p« + p' j p = + p* 

u=e(U +u*)'*'^v'+ew* 
X " y z 



When the perturbation quantities are introduced. 
Continuity ^ ^ ^^x ” ^ 



Momentum 




( 11 ) 

the resulting equations are 

( 12 ) 

(13) 
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Energy 




(14) 



We have formulated a 2D planar flow so that v’ = 0. To avoid writing 
subscripts, we now drop it in subsequent equations. The notation with 
subscript x means partial differentiation with respect to x* 

Starting with Equations (8) to (lO), one could obtain the linearized 
equations appropriate for heat and mass addition in axisymmetric flow. The 
equations have been solved for mass addition (sources) but not for heat 
addition in- axisymmetric flow. 

The solution for a source in two-dimensional planar flow will now be 
obtained. Consider a row of point sources evenly distributed along the 
z-axis. Axis orientation relative to the flow is shown in Figure 11. For 
a point P(x - x^, y - y^), only the sources within the Mach forecone influence 
point P as shown in Figure 11, To find the potential in two-dimensional flow, 
add up the sources along the z-axis from k to d. Points k and d are 
illustrated in Figure 11. From the geometry of Figure 11 




(15) 



Integration of a row of point sources yields the potential 



0(x,y,O) = - 




d 



(16) 



This integrates to 
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MACH FOREGONE 




FIGURE II. POTENTIAL FOR A LINE SOURCE IN SUPERSONIC FLOW. 



-1 






0(x,z) = - — sin 

Y(x - - p^(y - y^)' 

Inserting the limits gives 

0(5f,y) = ^ 

The potential is a constant. Think of a series of Mach cones with apex 
at the z-axis. Upstream of these cones, the potential is zero. Downstream 
of the envelope of Mach cones, it is constant. One can represent this by 

0(x,y) = - I |^(x - x^) + f(y - y^)j 

where I is the unit function. The velocity components are, for a unit 
source strength, 

and 

" If " ■ ^1^] 

The (+) sign in front of p gives right-running waves, whereas the (-) sign 
gives left-mnning waves. 



(17) 



(18) 



(19) 



( 20 ) 



( 21 ) 



The solutions for a heat source in supersonic flow are 

p. . (r 

pi = il “ l)Hq _ Py) _ lX-Z_Lla 5(y)l(x) 

2a p a U 



( 10 ) 



( 22 ) 

( 21 ) 

(24) 

(25) 
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One can verify that there are the appropriate solutions when m = 0 and Q / 0 
by substitution into Equations (12) to (14). The second term in Equation (25) 
is a wake. Pressure in the wake is the ambient value, and the velocity 
returns to freestream value. The solutions in Equations (22) to (25) are 
valid for all x and y including the source point x = y = 0. The symbol q, 
which has units heat/(length)(time), will be discussed shortly. The line heat 
source in a planar flow is illustrated in Figure 12. 

Similar solutions can be written for the case of mass addition; these are: 



= . ^ 6(x - py) 


(26) 


= ^6(x - py) 


(2'7) 


= 5(x - Py) 


(28) 


2 




= ^6(x - fy) 


(29) 



Substitution of Equations (26) to (29) into Equations (12) to (14) 
indicate these are solutions except at the mass source point x = y = 0. 

With Equations (22) through (29), one has solution for volume distribution 
of heat and in regions where m = 0. This is adequate to describe planar 
bodies with adjacent heat addition zones. 

Volume Heat Addition 

Since the heat may be added throughout a volume rather than a line source 

(the z-axis in Figure 11 ), it is worthwhile to formulate heat as a volume 

release. Define, as in Equations (lO) and (14) 

Q _ heat 

^ ( span ) ( area ) ( sec ) 
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FIGURE 12. SMALL PERTURBATION HEAT SOURCE IN A SUPERSONIC 
FLOW. 
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Then we define 



heat 

(span) (sec ) 
q = Qdxdy 

Another area element may be better than dxdy. Some element with its sides 
parallel to a Mach line is preferred. We think of Q = Q(x, y). Figure 13 
is a series of drawings illustrating the variables. 

To get heat per unit span 



L r-t 



H = 



0 



Qdxdy = QLt 

L 



(^ 0 ) 



if Q = constant 



Heat released in the complete volume is U = Hs. 

Tt * heat ^ 

sec (span) (time) 

Now to tie into the 2D flow equations, we use Equation (22) as an example, 

6(x - ^y) (31) 

At a point on a characteristic, u' results from a source on that same 
characteristic. If there is more than one source on that characteristic, 
one adds the influence of each. 



u' 



(y - l)(q3^ + Q2) 

2ypP 



6(x - Py) 



For many heat sources, in fact n heat sources, 

i = n^ 

- Py) 

i = 1 



(32) 



(33) 
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L 



♦ X 



(B) 





FIGURE 13. 

I LLUSTRATION OF QUANTITIES RELATED TO SPECIFICATION OF HEAT RELEASE. 
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The equation for u' is correct if one thinks of 



q. = Qdy = 



heat 



‘i ’ ' (span) (length) (sec) 
when the heat source is distributed in the volume. When heat sources are not 
discrete but are distributed, then one wants to integrate along a characteristic. 
Let s be distance along the characteristic. We express in terms of Q and ds 

q^ = Q sin p ds = dq (34) 

The perturbation velocity is 



s 



2YP^ ^ 



r s 



(y - 1) 



Q(x,y)ds 



0 0 
Since s,x,y are related quite simply 

X = s cos )Xf y = s sin x = fy 

one can integrate u' several ways: 

s y 

sin p f q|s cos ju, s sin jujds or j Q(py,y)dy 
■^0 0 



or 



Q(x,x/p)dx/^ 



(35) 



(36) 



Referring to the region boxmded by 0 < x < L and 0 < y < -t in Figure 13(b) 



and for case Q = constant inside the region 
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s y 

" r 

sin p Qds = Qt = j Qdy = Qt = 
The pressure perturbation would be 



X 



Q(x,x/f)^ = (tan ;u)PqI 



= Qt (37) 



0 



, _ (y - 1 )M0t 



P' = 



2a^ 



In terms of Q(x,y), the perturbation quantities are 



(38) 



u' = - 



(y - 1) 



Q(x,y) sin p ds 



(39) 



s 

' = Q(x,y) sin p ds 



(40) 



P' = 



(r - 

2af 



Q(x,y) sin p ds 



(41) 



P ' = 



_ (y - 1)M 



2a^p 



Y •• 1 

Q(x,y) sin u ds - — 

a U 



Q(x',y')6(y - y')X 



X* y* 



l(x - X ' )dx*dy' 

The wake is given by the integral above; as an example for x > + L 



(42) 



x' = X. + L 
1 



p' (y) = _ r_Z_l 
a u 



Q(x',y)dx' 



(4^) 



where the geometry is sketched in Figure 14. 
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FIGURE 14. GEOMETRY OF WAKE DOWNSTREAM OF A HEAT ADOPTION ZONE. 
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For a point inside the heat release region 






I 

Q(x',y)dx' 



a u 



(u) 



x' = 



The symbol x^ is defined in Figure 14. 

Heat Add ition at a Plane ; A Heat Front 

Consider a distribution of line sources which form a plane at angle 9 
to the flow, as illustrated by Figure 15. Terminology similar to that used 
in supersonic wing theory, i. e., subsonic and supersonic "leading edge," is 
employed to identify two cases, 9 ^ or ^ 1, The normal Mach nvunber 
is M sin 9. Since sin p = l/M, the normal Mach number is 

(«) 

You can see from an examination of Figure 15 that the heat source can influence 

flow upstream if M <T 1. However, if M > 1, the heat source does NOT influence 

flow upstream contrary to statements made by Tsien and Beilock.^^^^ Also it 

is apparent that point P is influenced only by right-running characteristics for 

M <1, whereas P is influenced by both families of characteristics for M > 1. 
n n 

The perturbation velocities for a heat front will now be obtained. 

Using Equation (22) and integrating along the heat front, the formula for u' is 

00 



u«(x v) = _ 

2ypP 



r 

5 j(x - 3 cos 9) - p(y - s sin 6^ 



ds 



(46) 



-00 



where s is distance along the heat front. Consequently, 
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( 47 ) 



^ (y - iJa zA 



2ypf p sin 0 - cos 0 



J 



CO 

s 

6[s - 

-<X) 






C03 6 “ p sin 9 



ds = 



(r - i)q 

2yp^(P sin 0 - cos 0) 



Since f = - 1, it is also true that f = l/tan fi. Using this fact 



I _ (y - 1 )q tan ju sin u 
2yp sin(p - ©) 



(48) 



For left-running waves, p. is positive; and for right-running waves, p is 
negative. For a subsonic "leading edge," 0 is less than p. Downstream of the 
heat source for < 1, Equation (48) indicates u^ will be negative; upstream 



of the source u^ will be positive. In addition 



(for M < 1) 
n 



“ii > i“ji 

where subscripts u and d indicate upstream eind downstream respectively. 

For a supersonic "leading edge," which is illustrated on the right-hand 

side of Figure 15, the left running Mach wave gives a u^ which is negative. 

The right-running Mach wave gives a u' which is also negative but which is less 

xl 

in absolute value than u£. Subscripts L and R denote left- and right-running 
Waves. Consequently, 

uj = <0 

Now let's look at v'; from Equation (23) one obtains for this case: 

- 00 

6 |^(x - s cos 0) - p(y - s sin 0)Jds 



v«(x y) = Iy. - 

^ 2yp 



(49) 



J 



-CO 
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Integration of the 6 -function yields 



]. = (y - -l )q sin M 

2yp P 6 cos 0 2yp sin (0 - p.) 

For the case < 1, is negative, and is also negative. 

Furthermore 



(50) 



V ' ■> V ' (for M < 1 ) 

I u I I d| n 

For the case M > 1, vJ is positive while vi is negative. In addition 

v^ = v^ - |v^| >0 (for M^> 1) 

which means the flow is turned toward the line source much the same as an 
oblique shock wave. For < 1, the flow turns away from the heat front. 

A natural question is, "Is there any special direction to the 
perturbation velocity vector?" Let v^ be the resultant perturbation velocity 
vector. It is given by 



V 

r 



-»• ^/ tan u sin u . tan u sin u\ , ^ I sin u sin u \ 

®x^\sin(0 - p) sin(0 + p) / ®y \sin(0 - p) sin(0 + p)J 



(51) 



The algebraic signs are correct for the case M^> 1 which means 0 > p. The 
symbol C equals (y - l)q/2yp. Let a be the angle of the perturbation 



velocity vector relative to freestream direction. Then 



sin u sin u 

tan a = sin(0 - p) si n(0 t 

tan u sin u ^ tan u sin aj 
sin(0 - p) sin(0 + p) 

Manipulation of Equation (52) leads to 



(52) 



, 1 sin u cos 0 1 

tan a = - t ^ — n = - 7 ^ 

tan p cos p sin 0 tan © 

The conclusion from Equation (53) is that the velocity vector is normal to 
the planar heat source. Since this is so, one can transform the supersonic 
case of heat addition to an equivalent ID problem. To do this, one uses 



(53) 
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du 



( 54 ) 



^lH 



(y. . 



u u 

where u = U sin 0. Proceeding by substitution of components from 
Equation (51) into Equation (54) gives 



du 1 

U sin 0 1 - c T 

n p 



(55) 



which is precisely the influence coefficient result when properly interpreted. 
See Table IV. For a subsonic ’^leading edge/* the velocity increases. This 
means that the flow is turned away from the planar heat source for < 1 
(analogous to Prandtl Meyer expansion) and toward the planar heat source for 
1 (analogous to oblique shock wave). 

As has been demonstrated, one of the ways of adding heat which lends 
itself to analysis is by means of a heat front. The example chosen by 
Tsein and Beilock^^^^ was a heat front. Oswatitsch^^^ ^ considers other 
aspects of heat fronts. 

Heat addition at a front is comparable in complexity to linearized heat 
addition. Combustion is a volume phenomenon. However, if the length of the 
heat release region is much less than a body length, a front is a reasonable 

(12) 

approximation. Detonation is a heat front of a very special nature. 
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A heat front with a subsonic leading edge located next to a wall will be 
preceded by an oblique shock. The oblique shock turns the flow away from the wall, 
and the heat front turns it back parallel to the wall. Recall that for a subsonic 
leading edge the heat front can influence the flow upstream. The corresponding changes 
in pressure are shown. The net change in pressure is positive although the heat 
front decreases pressure, 

A heat front with a supersonic leading edge will be followed by an expansion 
fan if the front ;Ls adjacent to a straight wall. Once again the net pressure 
change is positive. The heat front turns the flow away from the wall, and the 
fan redirects it so as to be parallel to the wall. 

Adding heat to a supersonic stream decreases the Kach number. When sufficient 
heat has been added, the normal component of the Mach number downstream of the 
heat front, M^, will become unity. The amount of heat required to make = 1 is 
known as ’’critical" amount of heat. With a heat front having a supersonic leading 
edge, heat can be added until =1, If more heat is added, the flow illustrated in 
Figurel6(b) changes to that shown in Figurel6(a). The oblique shock decreases 
the flow Mach number changing the leading edge of the heat front from supersonic to 
subsonic. 

Adding heat to a subsonic flow or a heat front with a subsonic leading edge 
modifies the flow upstream. One cannot specify the flow immediately ahead of the 
front. Critical heat addition does not have the same impact in the subsonic case 
since the flow upstream will be modified to adjust to the large heat addition, 

A detonation wave can be thought of as a shock wave followed by a combustion zone 
or heat front. In Figure. 16(a) if the oblique shock and the heat front had the same 
angle, it would be an oblique detonation wave provided the Mach number and amount of 
heat released were correct. 



40 








FIGURE 16 . HEAT FRONTS AND COMBUSTION FANS. 
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In the Prandtl-Meyer expansion, conditions are constant along a 
characteristic. Centered combustion fans are possible with conditions 
constant along a ray. However, the ray is not a Mach line, and the 
normal component of Mach number need not be unity. Combustion fans could 
occur in subsonic as well as supersonic flow. 

If M^, the normal component of Mach number, exceeds unity, the fan will 
turn the flow away from the wall like a compression. This is illustrated in 
Figure 16(c) on the left-hand side. If M^ is less than unity, the flow would 
be turned downward as shown in Figure 16(c) on the right-hand side. 

Starting with the continuity equation in cylindrical coordinates (r,f,z) 
and setting ®( )/8z = 0 and ®( )/9r = 0, one can derive the continuity 
equation for a combustion fan. It is 



dw^ 



dtnp 



1 dp 2 dp 



= 0 



(56) 



The velocity components are illustrated in Figure 17, Heat addition does not 
cause an irrotational flow to become rotational. One additional equation is 
the condition for irrotational flow which is 

dw, 

dp - "2 = ° (5'') 

The rate of heat release controls d-tnp/dp. One simple combustion fan occurs 
for d-tnp/dp = constant. For this case. Equations (56) and (5'^) can be 
readily integrated. 

There is a particularly interesting expression that can be derived for 
the pressure coefficient associated with a heat front in hypersonic flow. 
Consider a heat front at angle © relative to the main stream. The velocity 
components and angles are illustrated in Figure 18, The momentum equation is 
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FIGURE 17. GEOMETRY FOR COMBUSTION FAN. 
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FIGURE 18. DEFINITION OF ANGLES AND VELOCITY COMPONENTS. 
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(58) 



Pi - Pi^ = P2 P2^2 



which can be rearranged, since from continuity Pi^^i “ P2^2 

P2 - Pi ' Pl^l - Pa“2 - Pl^lS - P2> 



(59) 



From the geometry illustrated in Figure 18 

^2 ^ - ^2 
sinj^ - ©I sin 6 



(60) 



Noting that sin 6 and combining Equations (59) and (60) leads to 

- Pi = Pi'^i'^2 ® sin5 (61 ) 

If sin^0 1, heat addition has little influence on w or p. Hence 

^ For hypersonic flow Equation (61) becomes 



P 2 “ Pi = 6 sin 6 

and the pressure coefficient is 

c = 2 tan 0 sin6 (62) 

P 

This can be compared with the Newtonian expression which is 

0=2 sin ^6 (61) 

P 

Equations (62) and (63) provide sufficient tools to explore the external 
burning ramjet shown in Figure 19. Angle 5 is shown in the top figure and 
applies to Equation (63). 

Heat Addition and Vorticitv 

In this section it is shown that heat addition does not change an 
initially irrotational flow to a rotational flow. The curl of the velocity 
vector should be zero. For subsonic flow, one uses the equations from Tsien 
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-SHOCK WAVES 





FIGURE 19. EXTERNAL BURNING SCRAMJET. 
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and Beilock^^^^; and for supersonic flow, one uses Equations (22) and (23). 



-♦ _/aY' 8u* \ 

curl V = ej^ - 

Substituting Equations 22 and 23 into Equation ( 64 ) yields 

- fy) - P5')] * 



(64) 



curl V = e 



2np(x - fy) (x - fy) 



(65) 



Recall that - 1. In the above analysis, 5 has been represented by 



6(x - py) = expj^- n(x - py)^] (66) 

■ 

Since curl v vanishes, one concludes that heat addition does not generate 
vorticity. 

Mass Balance with Heat Addition 

To gain additional insight to diabatic flows, it is worthwhile to examine 
some simple flows using familiar methods. One such case is a mass balance. 

Consider a heat addition region as shown in Figure 20 with a control volume. 

Half of the control volume can be used for the mass balance. This half has 
width w. Writing a continuity equation for mass flux across control volume 
yields 

= (p * p’)v'(L - pi) + 2(p + p')v'pl + P^u^(w - 1/2) + (p^+ p^)u 1/2 (67) 

where p* and v' are averages of perturbation quantities along pi illustrated 

in Figure 20. The density pertxirbation in the wake is p'. Equation (67) 

w 

reduces to 

pv'L = - p^u^l/2 (68) 

when second order terms (e.g., p'v') are neglected. This can be rewritten as 
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) 




A A 
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V 



(69) 






P * ^ 

V 

p 2L 



Substitution of Equations (40) and (42) in (69) gives an equality, as one 
expects. 



EXTERNAL BURNING; THRUST MINUS DRAG 
Drag — Subsonic and Supersonic 

In subsonic flow of a perfect fluid around a body, there is no force on 
the body, either lift or drag. Motion of a body through a fluid would require 
work if there were a drag. This energy would appear in the fluid. For a 
perfect fluid, the perturbations in the fluid decay rapidly at large 
distances from the body. If the body has been moving for a long time and if 
there were drag, there should be significant motion of the fluid at large 
distances. 

In supersonic flow, there is significant motion at large distances even in 
a perfect fluid. The waves caused by body motion move out from the body eventually 
decaying to sound waves. Far from the body, the flow consists of outgoing sound 
waves. These waves carry away energy. Work overcoming drag supplies the 
energy. 

Relation Between Flux of Entropy and Drag 

(l3 ) 

It is well known (See, for example, K. Oswatitsch, *) that the flux of 
entropy through a control surface can be related to drag of a body within the 
control surface. To get a feeling for the extent of the waves causing drag, a 
sample problem was worked. The results are shown in Figure 21. 

Consider a 20° wedge at Mach 2, as shown in Figure 21. The wave geometry 
is shown and has been calculated using the method of finite waves. Width of 



49 




FIGURE 21. FLOW FIELD NEAR A WEDGE AT MACH 2 . 
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the base of the wedge is b. The lateral distance to a point on the shock 
wave is I* Distance is illustrated. The flux of entropy in and out of 
the control volume (neglecting the wake) is equal to entropy jump across 
the shock wave times the mass flux. The first expansion wave intersects the 
shock wave at point 2. The wave between points 1 and 2 accounts for 25*2 
per cent of wedge drag. The shock wave between points 1 and 3 accounts for 
39.3 per cent of the wedge pressure drag. It is necessary to go out 179 
base widths to obtain 99 per cent of the drag. 

The conclusion from this little study is that the waves are important to 
great distances from the body. In axisymmetric flow, the bow wave would decay 
more rapidly. 

Linearized Drag Due to Waves 

A wave is a mechanism whereby energy can be propagated. The energy takes 
two forms in the linearized or acoustic case. These are flow work and kinetic 
energy. It can be shown^^^ that an integral over space of the flow work equals 
the integral over space of kinetic energy. As an example, consider a half- 
diamond airfoil shown in Figure 22. The airfoil is moved exactly one chord 
length from position ABC to position DBA, and the undisturbed fluid is at rest. 
This motion does work Dc where D is the drag and c is the chord. The work appears 
as energy within the waves. Changes in the fluid motion occur normal to the 
waves. The disturbance originating at A now appears at A*, B at B’, C at C*, etc. 
The motion within the wave below A*B*C’ is due to earlier airfoil motion to the 
right of position ABC. The work done by drag moving one chord length appears 
within the volume, V, bounded by A*B*C*DKA. Consequently 
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FIGURE 22. DRAG EVALUATION USING WAVE ENERGY. 
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f 


Dc = 


l^p (U + u)^ + 


«> 


V 



( 70 ) 



where u and v are perturbation velocities. When there is heat addition, as 
illustrated in the lover part of Figure 22, Equation (70) changes. The force 
on the body is now thrust minus drag. Heat has been added. There is a 
left-running wave, AG, which moves into the fluid above the airfoil. Energy 
is added to the fluid bounded by ACG. Heat is added in a volume bounded by BFC. 
Downstream of this region there is a wake. Note that in the region EAC*B* the 
mass sinks cause u’ and v* to be positive* whereas the heat sources cause u'* 
and v'* to be negative. In region EAC*B* the mass sinks decrease p whereas heat 
increases p. The energy equation associated with the waves now is 



(T - D)c + (heat added) = 



I J* p[(U + U' + + (v* ♦ v")^] 



dV 



V 



C^l) 







f r 




Pd(^) ^ 


( energy )dV + 

. \J . 



V 



WAKE 



( energy )dV 
AGO 



Equation (71) is more conceptual than quantitative. There are some questions 
yet to be answered. Does the integral of flow work equal the integral of 
kinetic energy for waves caused by heat addition? Incorporating heat addition 
into the waves increases complexity. Can this method of drag and thrust 
evaluation lead to useful results and new insight? These questions will be 



*We use the notation that mass sources cause perturbation velocities u’ and v’; 
heat sources cause perturbation velocities u'* and v”. 
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considered in the section Wave Energy Due to Heat Addition which follows 
shortly. 

Linearized Momentum Relations 

In two dimensions, x and z, the drag formula in three dimensions is 
simplified. It becomes 



D = - 



j*(p - p^)dz - 
C 




(u 



u)[(U 



u)d 




The corresponding lift equation is 



( 72 ) 



L = 

where C denotes a closed contour surrounding the vehicle or wing. For 
small perturbations due to mass sources and sinks, the above equations 
become 



(p - 



P )dx 
^co 



pwjwdx - (U + u)dzj 



( 73 ) 



D = p 



00 



uwdx ~i 



2^co 



L2 2 

\J u 



w^jdz 



( 74 ) 



2 2 

Note that drag is due to second order terras, i.e. , uw, u , and w . First 
order terms drop out. For lift the small perturbation form is 



L = 






(udx + wdz) = p^U r 



( 75 ) 



In a following section the equivalent formulas for L and D will be derived 
for combined mass and heat addition. First, consider drag or thrust on a heat 
addition zone. 
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To illustrate the application of the heat addition equations, let us 
calculate the thrust of a line source of length, oriented parallel and 
normal to the stream. Equation (^’6) is taken from Chapter 7 of Liepmann 
and Roshko;^^^^ 



D = - 






pn^dA 



( 76 ) 



In their notation = x and = y. The control volume is shown in Figure 23. 
If Equation ('^6) is negative, one gets a thrust. Using the numbers for 
surfaces in control volume: 



D = - 



(p^U)(- U)dA 



-1 



- 






f (U - u')(p„ - P'„VE><- 



(H - 



u 



)(Po.- 



^ \'AVE^^* l)(v' )dA 



The integral over A^ can be integrated. The integral over A-, can be split 
into two parts, and the part with can be integrated. Equations (22) through 
( 25 ) were derived dropping terms in second order in v'/U, etc. The integrands 
of A^ and A^ reduce to Pa,UV. Consequently, Equation becomes: 
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FIGURE 23. CONTROL VOLUME FOR CALCULATING THRUST OF A LINE HEAT SOURCE IN A SUPERSONIC 
FLOW. 



0 = - P«,U^A^ + + 






J 



^P'VAKS^^ - 2 I P«,UVdA 



(78) 



We now need expressions for P'y^y^g which we obtain from 

Equations (25) and (23) respectively. Proceeding with p' 

I 






■ ; 6(y - y')l(x - x')ds 

a U 



(79) 



0 



where s is distance along the source; x and y, the coordinates where you want 
to evaluate Py^gJ and x* and y*, the coordinates of an element, ds, of heat 
source. When x> x*, 1=1; since we are evaluating A^-integral, x will be 
greater than x*. Also ds = dx* for the case at hand, 

.1 



- r) 

a U 



dx 



'■ <r-i)qt 6(yy.) 



a^U 



(80) 



Substituting Equation (80) into A^-integral of Equation (78), we get 



P'WAKS“ 



a U a 






(81) 



In a similar manner, an equation can be obtained for v* from Equation (23); 

t 



v'(x,y) = 



5|^(x - x') - ^(y - y*)]ds (82) 



Substituting Equation (82) in A^-integral of Equation (78) yields 
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(83) 






j 

*5 “ 



6[(x - x') - p(y - y')]dx*dx 



The coordinate system is chosen so that the heat source starts at x' = 0 and 
ends at x' = t. Also y' = 0. The coordinate of point c in Figure 23 is 
(^y^,y^), and the coordinate of point d is ( Py^ + '^•jy^). Consequently, the 
limits of integration are as shown 



py^ I i 

j J" 5[(x - X') - py 5 ]dx'dx 

py5 0 



(84) 



A change of variable helps evaluate the integral; let z = x - 



u( y - l)d 



I I 
r 



2yp 



6(z - x*)dzdx' = 



0 0 



(85) 



The integral in Equation (19) is in the standeurd form as shown in the 
Appendix so that it is readily apparent that it has value -t. Now combine 
Equations (78), (8l), and (85) to give 



D _ (y - l)qtU 
^ 2 
a 




2yp ) 



( 86 ) 



which vanishes since a = yp/p* Supersonic heat addition in an unconstrained 
(no walls) flow gives zero thrust or drag* 

Now look at the case where the heat source is of length t but normal to 
the streamlines. This is illustrated in Figure 24. Equations (78), (79), and 
(82) apply; however, ds is different as Is the limit of integration. Put 



the origin of coordinates at the center of the heat source. Equation (79) becomes 
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THRUST OPALINE SOURCE NORMAL TO FLOW. 






1/2 

{ S 
U J 



6 (y - y' )dy' 



(87) 



- 1/2 



The A^-integral in Equation (78) becomes 

r 

I y2 ix - l)q 1 5(y _ yt)dy'dy = (XJL-iMu 

J a U J a 

A, - 1/2 

since the integrations over y' and y overlap over the region - t/2 to 1/2, 
For the perturbation velocity, the equation for v' becomes 

M 



( 88 ) 



v'(x,y) = 



- ^(y^ - y'^dy' 



(89) 



^ 1/2 



where y^ is the coordinate for plane A^. The A^-integral is 

1/2 



r 



p^Uv'dy = 



J 



- P(y 5 - y')]dy'dx 



(90) 



A5 - V2 



The limits of integration for x are p(y^ - 't/2) to p(y^ ♦ t/2). To obtain a 
standard form for the integrand, we must eliminate ^ in the 6 -fxmction. This 
we do by letting z * My^ “ y* ) and dz = - ^dy'. Equation (90) is now 



u(:c--iJ^ 



® 2ypP 



5 (x - z)dzdx = 



P(y5 * 



(91) 
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The limits of integration for both z and x in Equation (91 ) are Hy. - t/2) 
to + ■t/2). Combining Equations ('^8), (88), and (91), one finds 



which is what we expected. 

It is interesting to draw the streamlines in the vicinity of the heat 
sources; see Figure 25. For the parallel heat source there is a conceptual 
difficulty which has a resolution. The "stagnation” streamline* C goes 
straight back since v' is an odd function and must be zero at y = 0. One 
can think of a streamtube bounded by streamlines B and D. Continuing, one 
can shrink this streauntube so that there is very little mass flow through it. 

Yet the streamtube diverges at angle 0. A finite amount of heat is added to an 
infinitesimal stream yielding a finite angle 0. Strange.* One avoids this 
difficulty by using volume heat addition. The case of the heat source normal 
to the flow is more easily digested. Streamlines F and J are acted upon by 
characteristics of a single family. Streamline G passes through the heat 
source into a region where it is influenced by both left-running (gives v’) 
and right-running (gives - v*) characteristics. In this region the streamline 
is not deviated from its original direction. Streamline G then moves into a 
region where only left-running characteristics act on it. It is deflected by an 
amount 0. Finally G moves into the wake where it accelerates and assumes its 
original direction. The speed along G is (U, U - 2u', U - u', U) as it moves 
from region to region. Streamline H is in the plane of symmetry; as a result, 
it does not deviate from freestream direction. 

*A body will stagnate the flow; however a heat source does not do so. 




( 92 ) 



61 



B 

c- 

D 



WAKE 





FIGURE 25. STREAMLINES IN VICINITY OF PARALLEL 
AND NORMAL LINE HEAT SOURCES. 
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Wave Energy Due to Heat Addition 

The model for determining wave energy due to heat addition is shown 
in Figure 26. A heat addition zone KFGH is translated its own Length L in 
time L/u^. Its new position is AFEL. The heat released is (hlL)(L/u^). 

Since there is no force on the heat addition region, no work is done moving it 
from KFGH to AFEL. The question of forces on heat addition zones was 
discussed in the preceding section. All of the heat energy released must 
appear in the waves and in the wake. 

In the time L/u^,, the disturbance at J has moved to point K. The 
disturbance at F has moved to B. The dashed lines AC and BD would be stream- 
lines if KFGH were stationary and the fluid moved by. The energy content of 
space boxinded by ANGEL is the same as that inside FJ1‘:HE. Part of the energy 
added appears in the region ABDC. This is the wave energy. The remainder 
appears in the wake in the region KFGH. From B to P, work is done on the 
fluid to compress it. From K to D, work is done by the fluid as it expands. 
The flow work cancels. We are now ready to write the energy equation. 



hll? _ 2 


1 kinetic energy] 


1 energy in \ 


u 


\ in ABDC / 


^ wake EFGH 1 



(93) 



The factor 2 comes from the fact that waves go both up and down in Figure 26» 
The kinetic energy in ABDC is 



2^co 



r» 




^ r 




2, , 1 














n i 



L L cos ^ 



u^ sin ju djdn 



(94) 



where dx and dy have been replaced by a coordinate system x = ^ and y = sin p. 
The symbol u^ represents the perturbation velocity normal to the wave. It is 
equal to 
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FIGURE 26. WAVE ENERGY DUE TO HEAT ADDITION. 
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2 ,2 ^ ,2 
u = u' + v' 
n 



( 95 ) 



where u' and v' are given by Equations ('^9) and (40) respectively. The 
magnitude of v' has been plotted in the upper portion of Figure 20. Inserting 
Equation (95) into Equation (94) and recognizing the functional form of v' 

(and u') as shown in Figure 20, one arrives at 



2KE in ABDC 
Lp^sin cos p. 



= 2 



lift 



it Mit 



- (u;^ ^ v;2) 



dj 



(96) 






The value of v' has been indicated in Figure 20. It is the maximum 
m ^ 

perturbation velocity. Integrating Equation (96) yields 

I 2KB in ABDC \ ^ ,2 ^ ^,2^ | ^ ^ 2 ^ ^,2j( _ = 

I 3' m m ' m m '' ^ ' 



Lp sin cos }i 



(9-’) 



■(“'m * 



The energy in the wake is 



/energy to\ , - T„) - p o 

^GH 00 p ' w ^ p cc 



(98) 



\wake EFGH/ p ''' w <»" p » T^ 

The distance FG is the length L. According to small perturbation results, 

the wake has freestreara velocity and pressure but reduced density and increased 

temperature. Knowing these facts and using p = pRT, Equation (98) becomes 

/ energy in\ t s 

^wake EFGh) p® ^w 

The density perturbation in the ;^ake p^ can be obtained from Equation (44). 

Substituting Equations (39) and ( 40 ) into Equation (9'^), substituting 
Equation (44) into Equation (99)» and combining the results in Equation (93) 
yields 
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(100) 



P„L tf> 

— = — (L - ~) 
^ 3 






00 



(y - Dot 



i2 



2yp 
WAVE ENERGY 



Y - 1 



♦ c T^Ln(} - 

p ® _2 



a u, 



CO 



WAKE ENERGY 



From Equation (lOO) it is apparent that the energy appearing in the waves is 

2 2 

second order. The wave energy was based on u' + v' which is second order, 

mm 

2 

The wake energy reduces to qIl /u^,. If Equation (lOO) is to be an equation, 
the wave energy must be neglected. 

There are some important conclusions from this exercise. First, the 
wave energy is second order. Second, most of the energy due to heat addition 
appears in the wake. Ehergy radiated away by waves is small compared to 
energy deposited in the wake. 

If one examines the linearized momentum equations in two dimensions, as 
given in the section. Linearized Momentum Equations, one finds drag is related 
to an integral of second order terms. As shown in the discussion of Figure 21, 
drag can be related to the waves generated by motion of a body. A conclusion 
is that wave energy resulting from either drag or heat addition is second order 
in the perturbation quantities. 

Control Volume Approach to Thrust Minus Drag with Heat Addition 

As an example of a two-dimensional planar flow, examine Figure 27. A 
half diamond has a region of volume heat addition indicated as a shaded area. 

The dotted lines outline an image heat addition region above the shaded region. 
Part of a control volume boundary is shown as a dashed line. Mach lines are 
drawn from the body to the boundary of the control surface. Below the boiindary 
are plotted the perturbation quantities for that location on the boundary. At 
the very bottom are the pressure, source strength, u* and v*. These are constant 
to mid chord and then switch to values opposite in sign. 
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FIGURE 27. CONTROL VOLUME APPROACH TO DRAG AND THRUST. 



The perturbation quantities due to heat addition are plotted about the right 
hand side of mid-page. Heat addition increases pressure, slows the flow, and 
deflects it downward. If the two plots were to scale, the positive pressure 
increment due to heat addition would exceed the negative increment due to flow 
expansion (mass sinks). The body would have positive thrust minus drag. 

The magnitudes of the perturbations due to heat addition are proportional to 
the length of the characteristic (Mach line) passing through the heat addition 
regions both real and image. Note that p', p", u', u", v*, and v" are of 
opposite signs. 

Let's look at drag and lift contributions due to the portion of control 
volume shown. In Equation (72) dz is zero. Hence the drag contribution 

D = J" p(U + u' + u'')(w' + w")dx (101 ) 

Multiplying the terms yields 

D = J* (pUw* pUw" + pu'w' + pu'w" ♦ pu"w' *• pu"w")dx (102) 

12 3 4 5 6 

The following summary gives the contributions (+), 0, or (-). 

Term 1 2 3 4 5 6“ 

0 to c/2 (-) 0 (+) 0 0 0 

c/2 to c {*) (-) (+) (-) (-) (+) 

The integral of pUw' over a full chord gives zero; however, the integral of 
pUw" makes a large thrust contribution. Terms 3 and 6 give drag, whereas terms 
2, 4, and 5 give thrust. 



68 



Lift is obtained from Equation ('^3). Once again for the segment of 
control volume illustrated, dz is zero. The lift contribution becomes 



r 

L = (p ' + p" )dx + 
1 2 



p(w'^ + 2w'w” + u"^)dx 



(103) 



The following summary gives the contributions of the various terms. 



Term 1 


2 3 4 5 


0 to c/2 (+) 


0 

0 

4 

0 


c/2 to c (-) 


{*) (*) (-) (*) 



The pressure p' averages to zero, whereas p'* gives a lift. The cross term 
w’w” reduces lift. 

To optimize the heat addition, the airfoil shape must be expressed 
analytically. Heat addition also needs to be described by some function of 
X and z. Heat release need not necessarily be uniform. 

Thrust from Cycle Point -of -View 

The thrust due to a propulsion device is given by 

F = m(w - w^) «• Ag(pg - p^) (104) 

If one evaluates the equation at the Trefftzplane, the pressure term may or 
may not drop out. This is one item which Oswatitsch^^^ ^ discusses. In 
Equation (104), F is thrust, and m is mass flow rate of fluid which receives 
a velocity increase w - w^ over the freestream velocity w^,. 

Consider three flows as shown in Figure 28. Heat is added at stagnation 
conditions and at some Mach number other than zero. The pressures at the 
exhausts are all equal, i.e., P 4 " P 2 ~ P 2 * Mollier Chart which depicts 
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(a) ISENTROPIC WITH NO HEAT 
AIX)ITION;THIS IS A REFERENCE 
FLOW. 



(b) ISENTROPIC WITH HEAT 
ADDITION AT STAGNATION 




(c) ISENTROPIC TO STATION ADD 
HEAT 2 T03; ISENTROPIC 3 T04. 



FIGURE 28. THREE FLOWS WITH TWO DIFFERENT WAYS OF ADDING HEAT. 
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the three flows is shown in Figure 29. It is apparent from energy considerations 
that 

dq + + (h^ ■ ( 105 ) 

Continuing 



- wp) = wdw = dq - (h^ - hp 



2 



Now 



dw > 0 
w, > wi 

The slope of a p = constant curve is 



^s i 



if* dq > (h. - hp = dh 



= T 



Combining Equations (106) and (108) leads to 



( 106 ) 



(Kl’) 



(103) 



5h 

wdw = dq - dh = dq - -^j ds 



The change in entropy of the flow is given by 

ds = ^ 



(109) 



( 110 ) 



A subscript q is used to indicate this is the temperature at which heat is 
added. Equations (109) and (llO) yield 



wdw = dq(l - |-) (111) 

q 

If one compares flow (a) with flow (b) of Figure 28, then the temperature 
ratio in Equation (111) is T^/T^^. If one compares flow (a) with flow (c), the 
temperature ratio is have dw > 0, one must have T 
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FIGURE 29. MOLLIER DIAGRAM FOR HEAT ADDITION. 
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As a nvimerical example, consider flight at M = 0.9. Heat is added at 
M = 0 in flow (b) and at = 0.5 in flow (c). Enough heat is added to 
accelerate the flow so that = 0.6. The same amount of heat is added in 

flow (b). The results are that and = 1.091. 

Equation (ill) is not usually applicable to supersonic or hypersonic 
flow. The entropy change in the exhaust was attributed solely to heat 
addition. In supersonic flow, shock waves cause entropy changes. Further- 
more heat addition may alter the shock wave geometry. 

Thrust and Drag in Various Flow Regimes 

For subsonic flow, the pressure term in Equation (IO 4 ) is zero. For 
supersonic flow, the pressure term if evaluated in the Trefftzplane is not 
zero; however, it is extremely small. To see that the pressure is not zero, 
examine the flow illustrated in Figure 30. The entropy in the Trefftzplane 
is greater than ambient; ds > 0. If dp were zero, then pwdw equals - Tds 
so that dw 0. Entropy is related to p and T by Tds = c^dT - dp/ p = c^dT. 
Consequently dT > 0. From p = pRT, it follows dp/p < 0. The continuity 
equation pwA must be satisfied. As shown in Figure 30, the streamtube area A 
does not change significantly. The changes in both p and w cannot be negative 
and still satisfy continuity. If dp > 0, then dw continues to be less than 
zero and dT "> 0. Also dp/p > dT/T. Consequently dp > 0. 

To account for the drag, the area A must be quite large. (See Figure 21. ) 
The streamtube influenced by a propulsion system will be small as shown in 
Figure 30. It is an area of approximately A^ in Equation (IO 4 ). The pressure 
at the Trefftzplane is due mainly to drag and is spread over area A, albeit not 
viniformly. The pressure term in the thrust equation is applicable to area A^. 
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FIGURE 30. DRAG IN DIFFERENT FLIGHT REGIMES. 
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Since thrust equals drag, the pressure term will be of the order of A /a. It 

e 

can be neglected for subsonic, transonic, and supersonic flow. 

For slender bodies in hypersonic flow, dw/w is small. If the thickness 

ratio is T= ( thickness )/(bo(Jy length), then dw/w is of the order T^. Heat 

addition also has little influence on dw/w, a fact that can be verified by 

examining Figure 7.3 in Volume I of A. H. Shapiro's The Dynamics and Thermo - 

( 7 ) 

dynamics of Compressible Fluid Flow . In hypersonic flow, heat addition 

causes large changes in pressure eind static temperature. 

In hypersonic flow, the shock waves are close to the body. The drag is 

concentrated in a small area, A^; see Figure 30. Thrust producing devices 

would influence an area comparable to Ajj. Consequently, the pressure term in 

Equation (104) must be retained as a significant factor. 

One can calciilate forces on a body by suitable integration of pressure. 

Pressure in excess of ambient on forward -facing surfaces causes drag; define 

this as pressure drag, D^. Pressure less than ambient on rearward -facing 

surfaces causes drag; define this as suction drag, D^. The ratio is 

shown in Figure 31 as a function of Mach n\jmber. Suction drag is very 

important in subsonic, transonic, and supersonic flow. However, as M 

becomes large and enters the hypersonic region, becomes insignificant. 

In fact, Newtonian theory assumes ambient pressure p^ is zero with the result 

D = 0. 
s 

A propulsion device will increase pressure on rearward surfaces as shown 
in Figure 31. Heat addition must increase pressure so that the area under the 
curve labelled "with heat" equals the area under the curve. Dp. In supersonic 
flow, the thrust device simiiltaneously cancels a large D^ and balances Dp. 
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FIGURE 31 CONTRIBUTIONS TO DRAG AND THRUST. 
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In hypersonic flow, the thrust device does the same thing except that in 
hypersonic flow is negligible. Heat addition in hypersonic flow balances 
Dp but counteracts only an insignificant D^. This has implications con- 
cerning propulsion efficiency. 

One more point is that in hypersonic flow is the major drag con- 
tribution, and this can be related to the bow shock wave. In supersonic flow, 
is important also so that the waves originating at the rear end of the body 
or at the rearward -facing surfaces are of equal importance to the bow shock 
wave in calculating drag. 

Contrast in Thrust at Supersonic and Hypersonic Speeds 

Consider a body immersed in a flow as shown in Figure 30. The force on 
the body in the direction of the freestream is 




A 



If there is no energy addition and if A is the area of the large streamtube, 

F equals D. (it has been assumed that expansion waves have made pressure 
increments negligible along the streamline from a to b in Figure 30. ) If there 
is energy addition to the stream and if A is, once again, the Trefftzplane in 
Figure 30, Equation (ll2) yields thrust minus drag. If one evaluates 
Equation (112) using the small streamtube in Figure 30 with flow m, then F is 
approximately thrust. * 

*It is not, of course, possible to separate out the thrust accurately using 
Equation (112). If it were possible, how much more simple the propulsion 
integration problem would be. 
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The continuity equation is, neglecting fuel added, 



0 = 



(pu - p^u^)dydz 



The energy equation in integral form is 



(in) 



2^2 
U ^ V 



u 



0 = 



pu(h + 
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•) - P + T 



dydz 



(114) 



In Equation (ll4), q is the heat added per unit mass, and v is the component 
of velocity normal to the freestream. By combining the preceding equations 
in the following way 



BO 

(u^ X thrust) + (hj^ q - ~) x continuity - energy = u^F 



one obtains 



= 

DO 



pu(h^ - h q) ♦ ujp - p^) - pu 



( ^2 ^ 2 
(u - «■ V 



dydz (115) 



Substituting h = e ♦ p/p into Equation (115) yields 







- p(u - U„) - pu- 



(u - u_) + V 



■ dydz 



(116) 



There are three terms in the integrand of Equation (116). The term pu(e - e^) 
is the flux of energy thrown away in the jet. The term pu [q - (e - ej] is the 
energy converted to work or jet kinetic energy. The term p(u - u^) is the flow 
work, and the last term in the integrand is the flux of kinetic energy of the 
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jet relative to the stationary surroundings. The left-hand side is the rate 

at which the thrust does work in overcoming drag. 

For transonic and supersonic flow, all terms are important. In hypersonic 

flow, (u - v/u^ 1. Changes in velocity due either to shock waves or 

heat addition are small. However Ap/p^, AT/T^, or Ap/p^ are of the order of 

unity or larger. Based on these conmients, one can neglect the jet kinetic 

energy in hypersonic flow. In hypersonic flow, the term p(u - u^) varies as 

2 

1/M as will now be demonstrated. Rewriting the term as 



p(u - u^) = p u (1 ^)(— - l) 

2 

one can introduce the Mach number dependence. In hypersonic flow, Ap/p 'v/ K 

OO 

where K is the hypersonic similarity parameter TM. The velocity perturbation 
2 

varies as T as discussed previously. Substitution of these quantities into 
the flow work term yields 



p(u - U^) 



P u 

^CD GO 



(1 + K^)(K^) 



(117) 



Compared to the heat addition term in Equation (116), the term p(u - u^^) is 
small for M 1 and K = 0(l).* Equation (116) reduces to 



F = 




(118) 



which is valid for K » 1 and K = 0(l), In hypersonic flow, Equation (118) tells 
us that the thrust is simply due to the heat added less the amount of internal 
energy thrown away in the exhaust. In supersonic flow, it is more complex; see 
Equation (ll6). 



^Requiring K = 0(1) is not a severe restriction. 
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APPLICATIONS OF EXTERNAL BURNING 



In the preceding sections some of the analytical tools for understanding 
external burning have been discussed. It is worthwhile to look at some 
possible applications of external burning. Four applications will be examined 
briefly in the following section; these are forces on a planar airfoil due to 
heat addition, transonic boattail and base pressure alleviation, base pressure 
modification in both the planar and axisymmetric cases, and spin recovery of 
aircraft. 

Before discussing the applications, a model for pressure rise due to heat 
addition near a flat plat will be discussed. This mode] uses some of the 
results from EXTERNAL BURNING; BASIC EQUATIONS AND SOLUTIONS and from 
EXTERNAL BURNING; THRUST MINUS DRAG. The model is a refinement of a similar 
model by Billig.^^^^ 

Model for Heat Addition Adjacent to a Flat Plate 

Figure 32 illustrates the geometry. There are three zones. Zone I is 
the freestream which is supersonic. Zone II is the heat addition region 
separated from Zone I by a heat front. Zone II is bounded by the wall, heat 
front, and streamline c. Zone III is downstream of the oblique shock com- 
prising the region below streamline c. The heat addition turns the flow by 
an angle 9. The streamline deflection angles and static pressures in regions 

II and ni must be equal along streamline c. In the following development, 

2 

terms in 0 are dropped as being negligible. 

The deflection angle is approximately: 
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The continuity equation is 

‘ ‘ ' “2“^2 
Across the oblique shock there is a pressure increase equal to 



( 120 ) 



P2 - Pi 









( 121 ) 



P 



Equation (121 ) defines Define Q' as the heat addition per unit mass of 
air flowing in the streamtube defined by streamlines a and c. The energy 
equation is 

or as 



°p”^T2 *^2^ ^p'^Tl 



( 122 ) 



( 122 ) 



(15) 



Using Equation (7,14) from Liepmann and Roshko, one can derive the 
following momentum equation 

- PjUjyj * f ■ '’/2 * ^ 1^2 ■ '’ 2 ®* ' ° 

Combining equations (119), (120), and (l2l) yields 



( 121 ) 



P2Va<“i ■ “2^ * 'P-! - '’2>'j'2 * = 0 



(124) 



A solution exists for u^ = u^ and p^ = p^. Let*s examine the consequences 
of that solution. From the continuity equation and equation of state 



P2U2y2 P^u^y^ 



(125) 



Since u^ = u^ and p^ = p^, this reduces to 




(126) 
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The energy equation becomes 



or 



m 12^, ^ 12 

2^2 ^23 " ° p ^3 * 2^3 



c T. «• Q‘, = c T- 
p 2 23 p 3 



(127) 



since = U 2 * From Equations (126) and (127) 



^3 ~ ^2 ^23 



^2 ''-^7 



c T, 
P 



(128) 



In the notation of the influence coefficients of Table IV, Equation (128) 
can be rewritten as 



dA dh^ 

T " 

P 



(129) 



In view of the equality expressed by Equation (129), the influence 
coefficients give the result dp/p = 0 and du/u = 0. 

The heat added by the heat front must increase pressure in the 
streamtube a-c by an amount equal to pressure rise across the oblique shock. 
Using this fact 



dp 



12 



^%^1 



P f 

Combining Equations (119) and (128), the heat added in Zone II is 
The difference between ^ order of 0* Hence 



( 1 ^ 0 ) 



(m) 



(132) 
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Total heat added to turn the flow by an angle 9 is 



= «12 * '^21 = ‘ f > d ”) 

The heat addition turns the flow by an angle 9 and increases pressure 
by the amount given by Equation (I'^O). 

Forces on a Planar Airfoil in Supersonic Flight Due to Heat Addition 

A problem which has been studied by several investigators is external 

(17 ) 

burning near a two-dimensional wing. One such study is that of Mager. 

Mager uses the linearized heat addition formulas, Equations (22) to (25), to 
find the forces on an airfoil in supersonic flight with heat addition in 
adjacent streamtubes. 

Consider an area bounded by and x^ in the x-direction and by a 
surface and a parallel side at distance h^. See Figure ?3. Within this area, 
add heat Q. From Equation (41 ) the pressure at the surface is 



-= 21 



^ - 1 QM \ r 

, 2 J 



sin u ds = 



(y - l)QMh 

app 



(lU) 



Where does the factor 2 come from? There must be an image heat source above the 
surface; otherwise, there would be a v' with flow through the wall. Combining 
Equation (134 ) and the definition of H [h is defined in the paragraph 



preceding Equation (3l).j gives 



Pw 



(y - l)HhM 
h. (x^ - x^ )a^p 



(135) 



i' f 

In Equation (135) the quantity h varies from 0 at x = x^ to h^ at x = x^ ♦ h^ cot fi 
It remains constant to x = x^, and then h decreases linearly to zero at 
X = x^ + h^ cot fi. This is illustrated in Figure 34, 
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FIGURE 33. HEAT ADDITION ADJACENT TO A SURFACE IN SUPERSONIC 
FLOW. 
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FIGURE 34. PRESSURE DISTRIBUTION DUE TO HEAT ADDITION 
ADJACENT TO A SURFACE. 
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It is of interest to examine streamlines near the surface in the region 
of heat addition; see Figure 35. The streamline passes along ABCD. It is 
shown dashed with an exaggerated slope below. The streamline is straight 
from the time it enters the heat addition region until it arrives at A. 

This is so since aA = a ’A. From A to B the right-running characteristic grows 
in length relative to aA, Hence the streamline is curved from A to B, From 
B to C the difference in length b’B - bB is constant; hence the streamline is 
straight but at an angle to the mainstream. From C to D less and less of the 
right-running characteristic passes through the image heat addition region, 
and V* becomes smaller and smaller until finally the streamline is horizontal 
at D, From C to D the streamline is obviously curved. 

When the streamline is moved closer to the surface, as is the case in 
the lower part of Figure 35, the deflection is considerably less. It is 
obvious that the streamline on the surface will not be deflected since 
aA = a*A = bB = cG = c*C = d*E = dD for this case. 

Comparing Figures 25 and 35, one gains some insight to the quandry 
discussed in connection with Figure 25. When the heat addition is confined to 
a mathematical line, as in the top of Figure 25, one has difficulty because of 
streamline deflection near the source. However, when the heat is distributed 
as in Figure 35, there is no problem. 

Based on Figure 34 and Equation (l35), it is apparent that one has in hand 
the tools necessary to study the forces on bodies due to heat addition, at 
least in the linearized case. Refer to Mager*s paper for performance of 
an airfoil with heat addition. 
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FIGURE 35. STREAMLINES IN REGION OF HEAT ADDITION. 
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Alleviation of Boattail Drag for a Plug Nozzle 

In the Mach number range 0,8 ^ Mq ^1.2, which is the transonic regime, 

the nozzle gross thrust coefficient decreases significantly, A factor which 

makes this transonic dip all the more important is the corresponding increase 

in the transonic drag coefficient. Typical data for a plug nozzle are shown in 

(l8 ) 

Figure "^6, which is copied from the NASA Memorandum by Harrington.'' ' Part (a) 
is of interest here while part (b) is not. Two flow phenomena cause the dip 

in One is the change in plug thrust. Figure 36(d) shows the behavior 

of plug thrust. Another influence is boattail drag identified in Figure 36(c) 

as primary flap drag. Plug thrust decreases and boattail drag increases in 

the transonic region. 

Is there anything that can be done to correct the loss of the transonic 

(19 ) 

Rabone shows that by going to small plug angles, some of the loss of 

plug thrust can be recovered, although there is still a dip near M = 1.2. 

(I8) 

Harrington^ uses a translating shroud that has a pronounced influence on 
the boattail drag; however, the loss of a plug thrust partially remains. 

An alternate approach is to burn externally on the boattail. The direct 
thrust produced by high pressure on the boattail has very poor SFC. However, 
a flow interaction may occur which greatly increases thrust. This phenomenon 
has been termed ’’wave trapping” by Fuhs. Figure 37 illustrates some essential 
features of the interaction. 

With heat addition on the boattail, the turning angle, 0, of the primary 
nozzle flow is less. See the angle 0 in the region labelled (l) in Figure 3*7. 

At the slipstream between the external flow and the nozzle flow, the expansion 
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Free-stream Mach number, Mq 

(a) Nozzle gross thrust coefficient. (b) Secondary total -pressure recovery requirements. 

Figure 3^. - Comparison of nozzle characteristics. Full-length plug; corrected secondary weight flow ratio, (Wj/Wpi^Tj/Tp = 0. 04. 

18 

(Reproduced from D. E. Harrington. ) 



Nozzle gross thrust coefficient, Cjg 




Ratio of secondary-air to free-stream total pressure, Pj/Pq 
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waves are reflected as expansion waves without heat addition. These reflected 
expansion waves further reduce the pressure on the plug. The reflection is 
indicated at point (2) in Figure T^. With heat addition, the expansion w^ve 
is reflected as a compression wave increasing pressure on the plug. 

In region (*^) without heat addition the pressure is low, and the Kach 
number is large. The oblique shock (4) has a shallow angle. It is necessary 
to match slipstream pTessure and flow direction for both the external and 
nozzle flows. At the match condition, the pressure rise across the obi ique 
shock is small, giving large expansion of the nozzle primary flow. 

With heat addition there are several changes to the flow in region (l). 
Pressure is increased. The flow is deflected to larger values of ©. The Mach 
number is lowered in region (l) due to heat addition. Either a weak or a strong 
solution of the oblique shock occurs. When the strong shock solution occurs, 
the external flow downstream of the oblique shock is high subsonic. Expansion 
waves reflected off the plug surface are not transmitted into the external 
flow. Hence waves are trapped and are reflected as compression waves. The 
compression waves maintain a high plug pressure. Having looked at the inter- 
action qualitatively, let us now examine the results of a sample quantitative 
calculation. 

W^;V 0 Trapping ; Quantitative Example 

(is) 

The same geometry plug nozzle was selected as that tested by Harrington. 

The primary flap forming the boat tail was turned inward at an angle of l'^^. A 
frees tream Mach number of 1.31 was chosen for the calculations; one reason for 
selecting this value of Mach number was to keep the external flow supersonic 
throughout when there was no heat addition. A nozzle pressure ratio of 6 was 
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chosen. According to Figure 5, at maximuin thrust a higxher NPR would be more 
appropriate for = 1.31. NPR of 6 is characteristic of cruise. In order that 
standard gas tables could be used, the ratio of heat capacities, Y> was set 
equal to 1.40. 

To obtain the flow field, a planar, two-dimensional, finite, wave calcula- 
tion procedure outlined in Section 12.9 of Liepmann and Roshko^^^^ was used. 

Wien the zone of influence, as defined by Oswatitsch^^ in Section 3.25 of his 
book, spans a distance small compared to flow radius, the two-dimensional planar 
procedure can give quite accurate answers for local values of flow. Near the 
plug noazle throat, this procedure should give small errors; however, near the 
plug tip, sizeable error should be expected. In addition, the change in wave 
angle at the intersection of finite waves was neglected. 

Results of the calculation are shown in Figures 38 and 39. Figure 38 shows 
the wave geometry. The dashed line is the slipstream between the external and 
the nozzle flows. 

Consider the wave* separating regions (l) and (2). It is reflected 

as an expansion wave from the plug, becoming an wave. The wave reflects 

from the slipstream immediately downstream of region (4). The value of the 

reflection coefficient is + 0.3. This means the reflected wave is E with a 

r 

strength of 0.3 of that of the incident wave. At the next reflection from the 
slipstream, the coefficient is ^ 0.12, which means continued expansion. At the 
next reflection, the coefficient is - 0.08; there is an extremely weak com- 
pression wave reflected from the slipstream. As a consequence of the expansion, 

*E means expansion, and C means a compression wave. Subscript r is for a 
right running wave, whereas subscript I is for a left running wave. 
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FIGURE 30. WAVE GEOMETRY FOR FLOW ALONG PLUG NOZZLE WITHOUT HEAT 
ADDITION. 




Q-| Q.^ 
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FIGURE 39. PRESSURE DISTRIBUTION ON PLUG. 



the pressure along the plug drops below ambient as shown in Figure 39. Pressure 
below ambient on the plug results in a negative pliig thrust. Integration of the 
curve shown in Figure 39 yields a value of the ratio, plug thrust to ideal nozzle 
thrust, equal to - 0.097. Direct comparison of this calculation with Hari'ington' 3 
data is not possible since all his tests at M = 1.31 were run with a shroud 
extending partially over the plug. However, the value of plug thrust has the 
correct algebraic sign and the correct order of magnitude. In fact, it is 
better than the correct order of magnitude; it is probably within 50 per cent 
of being correct. 

Now let's add heat to the region near the boattail as shown in Figures 3*7 
and 40. The length of the boattail is L. The thickness of the heat addition 
region is i. The pressure increase due to heat addition is obtained from 
Equation (41). Evaluation of the integral yields 



t = ( y - l)Mh^ 

^ ap 

Equation (136) can be solved for the heat addition h, BTU/volume second, 
then multiplied by hi. 



(136) 

and 



_ a£p»L BTU 

(y - 1 )M (length) (sec ) 



( 1 ?' 7 ) 



Multiplication by 2vr, i.e., the circumference of the boattail, leads to the 
time rate at which heat is released. Let H be the heating value of the fuel 
and m^, the fuel flow rate. The heat release rate is given by 

^ (138) 

The specific fuel consumption, SFC, is given by 



. „ 2nra£Lp' 



SFC = 



3600fhj 

AD 



(139) 
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HEAT ADDITION REGION 




FIGURE 40. HEAT ADDITION REGION AT BOATTAIL. 
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when AD is the change in boattail drag. The change in boattail drag is 
AD = Lp' sin a; the angle a is defined in Figure 40. Combining Equations 
(l38) and (l'’9) along with the expression for AD gives 



M = 1.92, H = 20,000 BTU/lb^, and a = 17°, and dividing by 778 ft-lb^/BTU 
leads to SFC = 1.685 per hour. This is not particularly an exciting value for 
SFC; however, for the strong solution the heat addition causes wave trapping and 
significantly increases plug thrust. Note that SFC is not dependent on the 
amount of heat added. Less heat gives less thrust; more heat gives proper-- 
tionately more thrust. 

Now let*s look at the flow on the plug. To start the solution, the slip- 
stream angle, 0, must be determined. Figure 41 is a plot of the curves showing 
pressure in the external flow as a function of 9 and the pressure of the nozzle 
flow also as a function of 0. The ordinate is pressure behind the oblique 
shock (or downstream of the expansion fan for the nozzle flow) divided by 
freestream stagnation pressure. There are two possible solutions; the strong 
solution was chosen for this example. Values are 0 = ■♦■ 0.8^, subsonic external 
flow, and a shock wave angle of 77°, Knowing 0 at the boattail-slipstream 
junction permits one to start the solution using finite waves. 

The wave geometry along the plug is shown in Figure 42. Since the 
reflection coefficients at the slipstream have the value - 1.0, the flow 
becomes periodic. The pressure along the plug, which is shown in Figure 43, is 
everywhere larger than ambient pressure. Integration of the pressure gives 



(y - 1)MH sin 6 

Inserting numerical values, a = 1000 ft/sec. 



3600ag 



( 140 ) 
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FIGURE 41. SOLUTIONS FOR INITIAL SLIPSTREAM FLOW ANGLE. 



100 




101 



FIGURE 42. WAVE GEOMETRY FOR FLOW ALONG PLUG NOZZLE WITH HEAT ADDITION. 



FIGURE. 43. PRESSURE DISTRIBUTION ON PLUG. 
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- plug thrust ^ Q 22 ? 

ideal nozzle thrust 

If heat is added so that pVprpQ = .027 on the boattail, the value of 3FC 
becomes less than 0.1 per hour J In Equation (l39), AD would now include 
both the decrease in boattail drag and the increase in plug thrust. 

Wave Trapping ; Discussion of Numerical Example 

In the course of making a numerical experiment, there are numerous 
assumptions and decisions to make. Having completed this part of an analysis 
of the influence of heat addition, one has new perspective. The decrease in 
boattail drag and the corresponding SFC are fairly straightforward calcula- 
tions, and the value of SFC = 1.685 is probably fairly accurate. 

Considering the plug thrust increment due to heat addition, it is apparent 
that a crossroads occurs at the point of taking a weak or strong solution. 

See, once again. Figure 41. For this example, the strong solution was taken. 

The strong solution gives a subsonic external flow at least initially. 

That fact, of course, is the basis for wave trapping. Somewhere along the plug, 
the external flow may accelerate from subsonic flow to supersonic flow. Once 
the external flow becomes supersonic, the waves are no longer trapped in the 
nozzle flow. Determination of the subsonic flow and the change to supersonic 
flow is extremely difficult. 

If one had taken the weak solution instead of the strong solution, the 
pressure distribution along the plug, at least according to current estimates, 
would be qualitatively similar to Figure 39. 

Changing NPR raises the curve labelled '*N02ZLE FLOW’* in Figure 41. 

As the nozzle-flow curve is raised, the weak and strong solutions converge to a 
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single solution at the point of tangency. See Figure 44(a) and 44 (b), 

This single solution is a strong shock solution. Increasing NPR beyond the 
point of tangency raises the curve for ’tVOZZLS FLOW'* so that an intersection 
of the nozzle flow curve and external flow curve does not occur. The oblique 
shock probably moves upstream to the point where the boattail is formed. 

This case is illustrated in Figure 44(c). This would give high pressure on 
the boattail but probably would give 3ow flug thrust. 

Adding heat on the boattail has two influences: (l) Due to change in 0, 

the curve for external flow in Figure 41 is shifted to the right; and (2) due 
to increased pressure ahead of the shock wave, the curve is moved vertically. 
This is illustrated in Figure 45 (a). For large NPR and small heat addition, 
there may be a solution as shown in Figure 45 (c). Additional heat release 
causes a shift in the curve making possible both weak and strong solutions. 

See Figure 45 (d). 

The validity of the flow shown in Figures 42 and 43 is doubtful due to 
the fact that the external flow probably accelerates from subsonic to super- 
sonic. Tne strong solution may not occur. Consequently, one is skeptical 
of the SFC = 0 . 1 . Additional analysis is necessary to obtain the plug 
pressure distribution with heat addition. Both the weak solution and refine- 
ment of the strong solution should be considered. One would expect ina heat 
addition to have a favorable influence on plug thrust. 

Base Flow Problem vi th Heat Addition ; Planar Flow 

One application of external burning which has been investigated both 
experimentally and theoretically is the modification of base pressure by heat 
addition. The heat may be added in the base recirculation zone, in the viscous 
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PRESSURE PRESSURE PRESSURE 




(a) TWO INTERSECTIONS; WEAK AND STRONG SOLUTIONS. 




(STRONG SAME AS ILLUSTRATED ABOVE) 



(b) TANGENT, STRONG SOLUTION. 






-OBLIQUE SHOCK 
^RECIRCULATION ZONE 





DISTANCE ALONG PLUG 



(C) NO INTERSECTION; SHOCK UPSTREAM OFBOATTAIL. 

FIGURE 44. 

INFLUENCE OF CHANGING NOZZLE PRESSURE RATIO WITH FIXED HEAT ADDITION. 
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9 

(a) INFLUENCE OF INCREASING HEAT ADDITION ON EXTERNAL FLOW. 




PRESSURE DISTRIBUTION 
AND FLOW SIMILAR TO 
FIGURE 44(c) 



(b) NO INTERSECTION WITH SMALL HEAT ADDITION. 
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PRESSURE DISTRIBUTION 
AND FLOW SIMILAR TO 
FIGURE 44(b) 



(c) TANGENT STRONG SOLUTION. 




PRESSURE DISTRIBUTION 
AND FLOW SIMILAR TO 
FIGURE 44(a) 



(d) TWO SOLUTIONS; WEAK AND STRONG. 



FIGURE 45. 

INFLUENCE OF CHANGING HEAT ADDITION WITH FIXED NOZZLE PRESSURE RATIO. 
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shear layer, or in the inviscid flow adjacent to the shear layer. In 1967, 

Billig prepared an excellent summary of work to date. More recently, 

(20) 

Roberts has written a survey on the subject. 

It is worthwhile to exaniine the base flow from an integral point of view. 
This is illustrated in Figure 46 . The control voliime for the momentum theorem 
is shown as dashed lines with various surfaces identified by a, b, c, . . . . 
Using notation of Chapter 7 of Liepmann and Roshko,^^^^ the momentum theorem is 



pu, ( . 



i.n. )dA 
J J 



jpn, 



dA = 0 



(141) 



One can add an integral of ambient pressure over the control volume. For 
base pressure i = 1. 






(p - p^)n^dA = 0 



(142) 



For the control volume shown in Figure 46, there are four integrals. 



(p. - 



p^)(- l)dA + 



ab 



pu^(- l)dA + (p - p^)n^dA + 



piL dA 



= 0 (141) 



be 



cd 



de 



<Pb ■ p» 



ab 



de 



pu^dA + - pu^dA ♦ J (p - p^)n^dA (144) 

O'"©"© 

Heat addition in the inviscid flow above the edge of the shear layer will alter 
the shape of the streamline from c to d. Changing streamline cd will not 
influence integral above. It is not readily apparent how integral 
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EDGE OF SHEAR LAYER 



changes due to modification of streamline cd. Of course 
directly influenced. 

To gain some insight into the anticipated change in 
can use the continuity equation which is 



integral 




is 



integral 



one 



pu^dA = 



pu^dA 



(145) 



be de 

Since integral over surface be ioesn*t change when streamline cd is changed, 

this means the integral over de does not change. Any variation in integral 

2 

of Equation (144) is due to variations in p and u^; however, Equation (145) 
tells us any change in p is balanced by a compensating change in u^. 

If the flow above the edge of the shear layer illustrated in Figure 46 
is inviscid and if only one famdly of characteristics is significant, then 
one can relate local streamline slope to pressure. To demonstrate the influence 
of the shape of streamline cd on base pressure, a sample problem was worked. 

This is illustrated by Figure 47 . Streamline cd is used with the same lateral 
displacement but with large turning angle in the upper diagram. The pressure 
distribution along the streamline is shown in the graphs below each streamline. 
From the curves one can evaluate the integral or equivalent series 



I = 




n,dA = ) IJ-AA. 

1 / Pt J 



( 146 ) 



where AA. is the element of area projected in the direction. Numerical 
values for I are 12.8 for upper curve and 21.7 for lower curve. Relating 
Equation (146) to Equation (144) gives 
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FIGURE 47. PRESSURE ON EDGE OF SHEAR LAYER. 



no 






(147) 



Prp 




(p - p^)n^dA = p^(l - 52.8) 



cd cd 

The sharp turning curve gives a value of - 40p^ for integral of 
Equation (144), whereas the gradual turning curve gives - 31.1p^, The base 
pressure woul.d be greater with the gradual turning streamline cd* It is 
obvious from this analysis that one wants to maintain a high pressure on 
streamline cd to increase base pressure. The heat addition zone shown in 
Figure 48 will do that. 

Now let’s return to the problem of how much pressure increase results 
from the heat addition zone. Consider a heat addition zone located above a 
slipstream as shown in Figure 49. This is a model to represent the flow 
depicted in Figure 48. One replaces the lateral gradient in Kach numfoer 
with a slipstreami having supersonic flow above and subsonic flow below. V-ith 
a lateral Mach nmber gradient, as in Figure 48, the characteristics reflect 
from the sonic line in the form of a cusp. 

A paragraph will now be used to explain the notation E^, C^, E^, 

The syiTibols E and C are for Expansion and Compression waves. The subscripts 
r and I represent right running and left running waves respectively. 

In Figure 49 the heat addition zone has a leading edge swept at angle |i. 
This gives a finite wave C^. The pressure alorig the slipstream remains at 
pressure p^. The waves from the heat addition zone and the reflected waves 
turn the flow while keeping pressure constant. One way to describe this is to 
say that the compression waves from the heat addition zone have cancelled the 
expansion waves which would have originated due to kink in streamline at 
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HEAT ADDITION 
ZONE 




FIGURE 48. 

INTERACTION OF HEAT ADDITION ZONES WITH BASE FLOW. 
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FIGURE 49. INTERACTION OF A HEAT ADDITION ZONE WITH A SLIPSTREAM. 



point (^. The turning angle is 2© at point or twice the turning angle due 
to C^. The reflection has effectively doubled the influence of the waves from 
the heat addition zone. 

The doubling of effectiveness can be seen in the following analysis: In 

free unbounded space the turning angle is given by 

V ' _ (y - 1 )hb 



tan 9=0 = 



(148) 



U ♦ u' U 2ypU 

Equation (148) is obtained from geometry and Equation (40). From Equation (41) 



, _ (y - l)Mhb 



P' = 



2a p 



(149) 



where symbol b is defined by Figure 49. Due to turning of flow by an 
amount 29, the increment of pressure decrease would be 



p" = 



_ 2yM^p9 



( 150 ) 



From the flow field analysis illustrated in Figure 49, the following 
equation must be true 

2p' ♦ p'* = 0 (151) 

Combining Equations (148), (149), (150), and (l5l), one can verify that 
Equation (I5l) is in fact true. 

The iriechanism for altering base pressure by heat addition in the adjacent 

inviscid flow is apparent if one retraces the discussion associated with 

Figures 46, 47, and 48. To siuninarize, the streamline cd must turn inward 

toward the axis or plane of symmetry. Without heat addition, this turning 

causes the pressure to drop. Sharp turning gives low pressure and low base 

pressure; see Equation (144) and Figure 47 . Heat addition allows this inward 

( 21 ) 

turning without loss of pressure as shown in Figures 48 and 49. Strahle’s 
analysis suggests a new way to overcome base drag. 
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Base Floy Problem with Heat Addition ; Axisymmetric 




planar geometry, i.e., wing-like objects. Projectiles are, of course, 
axisynimetric. This section examines some of the changes to be expected when 
the filanar results are replaced by calculations for the axisymmetric case. 



flow field for axisymmetric heat addition. To verify the procedure for 
handling the characteristics near the axis, a model of a spherical, radially 
expanding flow was used. The results of that calculation are shown in Figure 50* 
Since the flow is known, one knows Kach number at any radius R. The l^fech 
number calculated from the characteristics solution is shown in the third 
column of the data summary in Figure 5Q» The Mach number based on the value 
of R and flow area is shown in the first column. 

A planar radial flew has an area variation linear in R, whereas the 

2 

axisymmetric radial flow has a R dependence for area. This can be represented 



Partial solutions of the flow field with heat addition were obtained for 
two Cases. Both cases had an annular heat addition zone as illustrated in the 
lower left-hand corner of Figure 51. The heat addition was sufficient to cause 
a maximum turning angle of 0 = - 6^ at a radius of 6 inches. In the first c; 

the heat zone was bounded by 6 ^ r ^7 and 0 ^ z ^ second case, it 

was 6 ^ r ^8 and 0 ^ z ^3. For the first case, shown in Figure 51, = V5. 




by 




(152) 



Planar Axisymmetric 
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FIGURE 50. CHARACTERISTICS SOLUTION FOR AXISYMMETRIC 
RADIAL EXPANSION, 
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ANNULAR HEAT ADDITION 




For the second case, was chosen as 2, and the radial extent of the heat 
addition zone was increased to spread out the waves. This was done to avoid 
merging of the characteristics to form a shock wave. 

Pressure versus radius is plotted in Figure 51- The mid-characteristic, 
which has points 1, 2, and 3, has an imperceptible increase in pressure as 
radius decreases. The characteristics converge to foi*ra a shock wave. When this 
happens, the techniques outlined in the book^"^^^ and NACA TN^^^^ by Ferri are 
apjdied. 

Along the characteristic, which starts at r = 6^’ with 9 = - 6^, the 

pressure rises slightly to a radius of r = 2'\ Between 0 ^ r ^ 2, the pressure 

increases rapidly. Due to the large mesh size, the pi*essure variation at 

points 12, 13, and 14 inay not be correct. Certainly the end point, number 14, 

is correct. The dashed line represents the curve that would probably 

obtained with finer mesh size and a more accurate treatment of the shock waves. 

2 

Note that for = 2, even a normal shock wave is nearly isentropic, which 
permits treatment of shock waves as finite isentropic compression waves. 

The conclusions from Figure 51 are as follows: 

(a) In the axisymmetric geometry, compression waves tend to merge more 
quickly. 

(b) The increase of pressure with decreasing radius is very small until 
the axis is approached. 

Now let’s look at the second example which was chosen to avoid a shock wave. 

The characteristics for the second case are shown in Figure 52* Increasing, 
the spacing of the waves did not completely avoid a shock wave. A normal shock 
extends from points 13 and 14 to the axis. 
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X — FIGURE 52. INITIAL CHARACTERISTICS IN REGION 

HEAT ADDITION ZONE t ' INTERNAL TO ANNULAR HEAT ADDITION, 



I 




I I Q 



The pressui*e as a function of radius is shown in Figure 53 for the mid-wave 
(wave starting at r = 6 with © = -3°) and the terminal wave (wave starting at 
r - 6 with 0 = - 6^). As before in the previous example, the pressure does not 
increase significantly until very close to the axis. This is time for botti the 
mid-wave and the terminal wave. The pressure is shown as a dash-dot-dash lii*e 
between point 14 and the axis since this was not calculated in detail. 

Figure 54 gives information about the streamlines and flow deflection. 

Right on the axis, of course, the streamline is the axis. There is large 
turning near the axis giving small radius of curvature for the stream.! jne. 

This balances the large radial pressure gradient near the axis. 

The conclusions from the second case are the same as the first case. 
Preliminary Coirments Bas ed on a Cui^sory Look at Axisyi^imetric Case 
The change in base drag is given by 

AD^ - (151) 

Each of the two terms, Ap, and A, , will now be discussed. The results shown 

b b 

in Figures 51, 52, and 53 suggest that Ap^ may not be increased significantly 

as a direct result of changing geometry from planar to axisymirietric. However, 

the results shown in Figures 51, 52, and 54 more closely apply to pressure at 

the dividing streamline and not to base pressure. The flow in tlie base region 

may magnify the pressure, just as the initial pressure due to heat addition at 

r = 6 is increased slightly due to axis;>Tnmetric geometry. 

Compression originating with the heat addition will be reflected at the 

sonic line. See Figure 48 and point of Figure 55. These compression waves 

will be reflected from the wake of the heat addition zone; see point of 

Figure 55. The reflection may be negative, i. e. , changing the compression waves 

to expansion waves. This negative reflection would not be favorable. 
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FIGURE 53. PRESSURE DISTRIBUTION ALONG CHARACTERISTICS OF FIGURE 52. 
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SMALL RADIUS OF 
CURVATURE BALANCES 
HIGH PRESSURE 



(A) STREAMLINES NEAR AXIS 



RADIUS 




(B) FLOW DEFLECTION AS A FUNCTION OF RADIUS 



FIGURE 54. STREAMLINES AND FLOW DEFLECTION. 
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FIGURE 55. BASE PRESSURE INCREASE DUE TO HEAT ADDITION. 



Now let’s look at the area terra. Consider the same amount of heat, Q, 
released in both planar and axisyrametric case. In the axisyraraetric case, Q is 
distributed along a circle of radius giving A/2^r^f BTU/length sec. In the 

planar case, Q is distributed along two lines of length ^^^2* line is at the 

top, and the other at the bottom. See Figure 56, For the planar case, the 
heat release is Q/^rr^, BTU/length sec, which is identical to the axisyrametric 
case. The area influenced by the heat release is in the ratio 

base area; axisyraraetric case “^r? r 

— i = - 1— (154) 

base area; planar case ^^1^2 ^^2 v a.; 

Making ^r^/r^ small helps increase pressure at r^ due to heat addition; see 
Figure 56. Making makes the base area very small in axisyraraetric 

case as compared to the planar case. One factor helps; the other hurts. 

As a final remark, it should be stated again that these are preliminary 
thoughts and conclusions. 

Spin Recovery Wing Heat Addition 

Some modern fighter and attack aircraft can be flown into spins of 
such nature that special equipment is needed for recovery. The drag 
parachute may be deployed to lift the aircraft tail. Before normal flight 

I 

is restored, the parachute must be jettisoned. An alternate method of 

I 

obtaining a torque on the aircraft would be heat addition. Figure 5^ illustrates 

I 

► 

a modern fighter in a spin. The engines run normally and produce thrust. By j 

1 

spraying fuel into the exhaust stream, the heat addition deflects the jet down- 
ward. A torque is developed which tends to lift the tail. Let's compare a , 

parachute recovery with heat addition. 



124 




FIGURE 56. COMPARISON OF BASE AREA FOR PLANAR AND AXISYMMETRIC GEOMETRY. 
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FIGURE 57. SPIN RECOVERY USING HEAT ADDITION. 
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Assume the aircraft has a moment of inertia of 6 x 10 lb ft about the 
pitch axis. Assiime the installed thrust is ?5,000 lb with afterburner. If 
the flow can be deflected by 6*^ with heat addition, the force tending to lift 
the tail is 4000 lb. If the aircraft has a moment arm of 22 feet, the 
restoring torque is 8.8 x lO^lb ft. Time required to point the aircraft 
downward, i.e., tui’n 90*^ in pitch, is 



t = 




46 sec 



If the sink rate is 140 ft/sec, the recovery altitude is 64 OO ft. 

One can estimate the amount of heat to be added for a deflection of 6° 
using Equations (119) through (129). The temperature of exhaust must be 
increased by 30 per cent. 

Now consider a parachute of 10 ft diameter with a drag coefficient of 2. 
The drag at a sink rate of 140 ft/sec will be 

D = = 4200 lb 



This is comparable to the force due to thrust deflection. Performance of 
parachute and heat addition is comparable. 



SUMMARY AND CONCLUDING DISCUSSION 

Use of an area rule or quasi area rule implies optimization. The area 
rule can be used to translate an optimum axisymmetric body, e.g. , von Karraan 
ogive, to a three-dimensional body. The optimum axisymmetric body shape is found 
by application of the calculus of variations to the source distribution function, 
fl^x). When there is heat addition, there is a new function to optimize, Q(x,r); 
see Figure 55. Both f(x) and Q(x,r) must be simultaneously optimized. The 
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axisyTTiinetric, small perturbation, solution for heat addition has never been 
obtained; it is an essential building block of the quasi area rule. 

Wlrien the body develops lift with aerodynamic circulation, there is another 
complexity to the analysis. Vorticity due to lift must be included in the analysis. 
There may be an optimum balance between lift generated aerodynamically and lift 
generated by heat addition. 

To af'ply the optimization teclinique to turbojets, ramjets, or other 
internally burning, air breathing, engines, a different analytical technique 
is required. A combined thrust-drag technique was briefly Introduced. This 
method uses an energy disc. 

By discussing both one -dimensional and two-dimensional, planar, heat 
addition, the connection between heat fronts, combustion fans, and one -dimensional 
heat addition was demonstrated. Heat addition in the planar case may be represented 
either by line sources or by a volume distribution function, Q(x,r). In several 
ways heat addition is similar to the flow generated by a solid body. The 
pressure cn the flat plate of Figure 32 is the same as that caused by a wedge 
of angle 9. Pressure and velocity perturbations are identical for both solid 
body and heat addition. Differences exist between flow over solid bodies 
and flow through a region of heat addition. One difference is drag. 

There are three approaches to calculate thrust minus drag in an inviscid 
flow: integration of pressure over the body surface, momentum control vol'.une, 

and wave drag. Integration of j)ressure is the most direct of the methods. 

Figures 22 and 27 illustrate the wave energy and momentum approaches for a 
thrusting, lifting, planar airfoil. Figures 20, 23, and 26 illustrate the wave 
and nioment'Jjm approaches for an isolated heat addition zone in an infinite medium. 

The merit of a particular method depends on the geometry. 
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Applications of external burning include lifting, thrusting airfoils, base 
or boattail pressure modification for an exhaust, development of base thrust 
for projectiles, and spin recovery. External burning for airfoils and projectiles 
has been studied extensively both theoretically and experimentally. A new 
phenomenon occurs for the exhaust problem; this is wave trapping. A 
preliminary analysis indicates heat addition for spin recovery is comparable 
to use of a parachute. 
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